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ABSTRACT 

We use hydrodynamic simulations of minor mergers of galaxies to investigate the nature of surface 
brightness excesses at large radii observed in some spiral galaxies: antitruncated stellar disks. We 
find that this process can produce the antitruncation via two competing effects: (1) merger-driven 
gas inflows that concentrate mass in the center of the primary galaxy and contract its inner density 
profile; and (2) angular momentum transferred outwards by the interaction, causing the outer disk to 
expand. In our experiments, this requires both a significant supply of gas in the primary disk, and 
that the encounter be prograde with moderate orbital angular momentum. The stellar surface mass 
density profiles of our remnants both qualitatively and quantitatively resemble the broken exponentials 
observed in local face-on spirals that display antitruncations. Moreover, the observed trend towards 
more frequent antitruncation relative to classical truncation in earlier Hubble types is consistent with 
a merger-driven scenario. 

Subject headings: galaxies: kinematics and dynamics, galaxies: interactions, galaxies: formation, 
galaxies: evolution, galaxies: spiral, galaxies: structure, methods: n-body simula- 
tions 



1. INTRODUCTION 



Th rough the pione e ring s tudies of IPattersonI (| 19401 ) 
and Ide VaucouleursI ([1959), it was first recognized 
that the stellar disks of most spiral galaxies are well- 
approximated by an exponen t ial su rface brightness pro- 
file. However, Ivan der KruitI (|1979( ) later found that in 
many cases this description breaks down at large radius, 
whe re the disk surface dens ity appears to be truncated 
(see lPohlen fc Truiilloll2006l and references therein). 

Many studies of truncated disks have examined 
the s urface br ightness profiles of local edge-on sys- 
tems (jvan der Kruit 1979; Barnabv & T hronso n 1992,; 
Barteldrees fc D ettmar 1994; Pohlc n et all I2OOOI: 



de Griis et aLri2001 : van der Kruit .20011: iKregel et all 
2002t iFlorido et al.l l2q06allbl) and at h igher redshifts 
z ~ 1 (iPerezI 120041 : iTruiillo fc PohlenI l2005h . This 



choice of inclination facilitates detection of truncations, 
but is subject to potential biases owing to the effects 
of d ust extinction and line -of-sight integration (see, 
e.g., iPohlen fc Trujillol l2006t ). Observations of face-on 
systems, which mitigate these complications, were also 
successful at detecting truncation, but found that the 
surface brightness profiles are better repr esented by a 
broken exponential than a hard break (iPohlen et all 
200I lErwin et all 120051: IPohlen fc Truiilld 120061) 
Hunter fc E lmegreenl (j2006f ) also note the existence of 



double-exponentials in their sample of very late-type 
spirals and dwarfs. 

These studies also uncovered a broad range of behav- 
iors, including some disks that fo llow a pure exponential 
profil e out to very large radiu s (|Bland-Hawthorn et al 
120051 IPohlen fc Truiillol l2006f) . and - as lErwin et al 
( 2005f ) first observed - some that are antitruncated, with 
an excess surface brightness relative to a n ex 
tial profile fitted to the inner disk ( Erwin et al 
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|2007() . These extended stellar disks dominate the light 
past 4-6 scale lengths, and have flatter profiles with scale 
lengths that are ^ 50% larger than the inner disk. 

A number of authors have proposed theoretical expla- 
nations for truncated disks. Dynamical arguments for 
star-formation thresholds have been s uccessful in mo- 
tivating the l ocations of truncations (jKennicuttl 119891 : 
iSchavj I200I iNaab fc OstrikeJ l200l). However ob- 
servations of ultraviolet emission Irhilker et all [200l 
iGil de Paz et al.l l2005. 2006; Boissicr et alJ|2006D, young 
stellar populations (Cuilland re et al., .200ll). and HII re- 
gions (jFerguson et a l. 1998) in extended disks suggest 
that there is indeed star formation occurring beyond im- 
plied thresholds. More recently, star formation models 
including a variety of triggering mechanisms (gravita- 
tional instabilities, spir al wave shocks, and stellar and 
turbulent compression; lElme green fc Huntei]l2006h and 
N-body simulations of angular momentum redistribution 
via bar instabilities (Debattista et al. 2006) have been 
more successful at explaining classically truncated disks. 
Despite these successes, theoretical mechanisms for pro- 
ducing antitruncated disks have received comparatively 
little attention. 

While it is clear that secular processes can influence 
disk structure and may produce truncated stellar 
disks, it has also become apparent that disk galaxies 
exist within a hierarchical unive rse, in which merg- 
ers are a frequent oc currence (iLacev fc Cold Il99l 
ISomerville fc Kolattl[T99 9: Somc rville et al.ll2000D . Fur- 
thermore, these mergers are likely to play an important 
role in shaping the appearance of galaxies. This is 
certainly true of collisions between spiral galaxies of 
equal mass, so called major mergers, which have been 
suggested as the dominant formati on mechanism for 
prese nt -day el li ptical galaxies (iToomre fc Toomrd 
1971 iTtaomrg Il977l: iNeeroponte fc Whitd 119831: 



Barnes! 119921: Hernquist fc Weinberg 19921: Hernguistj 
1998hl : iSilkfcWvsd 119931 : INaab fc BurkerU l200a 



iRobertson eralll2006allbl : ICox et al.|[2006b[ ). 

There has also been considerable study of the 
effects of minor mergers {Mprim/Mgec ^ 3) 
on the ver tical structure and dyna mics of stel- 
lar disks (iQuinn fc GoodmanI IQ i: 



Sellwood et al 



Font et al 



200lt 



et al.1 



1991 iWalker et all llQQ fi t iHuang & Carlberd 11997 
' 19981: IV elazquez fc White 199^ 
Ardi et al. 2003; Brook ct al. 200l 



20051120061: iGauthier et JB OOG: Hayashi & Chiba 200 



Kazantzidis et al.ll2007f ). in addition to observational ev- 



idence for past interactions with satel lites as t he origin of 
the Milky Way's thick disk (iFreeman fc Bland -HawthornI 
120021 : iGilmore et al.ll200l IWvse et al..,2006, ). Moreover, 
tidal structures indicative of recent min or mergers have 
been observed in both t he Milky Waj^ (Newberg et al' 
20021: llbata et al " 



McConnachie et al 



200l and 
120031) . 



M31 pbata et al. .20011 



In this work, we explore the effects of minor merg- 
ers on the structure of stellar disks at large radius and 
find that, under certain conditions, this provides a viable 
physical mechanism for producing antitruncated disks. 
We demonstrate this process using a set of hydrodynam- 
ical simulations, which are described in §[51 An overview 
of the merger process and specifically the dynamical re- 
sponse of the stellar disk during a minor merger is pro- 
vided in § [31 § [31 summarizes the surface density profile 
fitting procedure, and §[5l and [SI describe the dependence 
on parameters of the interaction. In § [3 we discuss our 
results in comparison to observations of antitruncated 
disks and within the context of hierarchical galaxy for- 
mation. Finally, we conclude in § [SI 

2. THE SIMULATIONS 

For this study, we consider the effects of a 1:8 
merger on the stellar surface density of the primary 
component's disk. These intera ctions are both cos- 
mologically common (see , e.g.. [Lacev fc Cold 119931 : 
iSomerville fc Kolatti 119991: ISomerville et al.l I2000D and 
kinematically important enough to play a significant 
role in determining the appearance of most present-day 
stellar disks, while at the same time largely preserving 
the overall disk structure (see, e.g.. iQuinn et al.|[l993l: 
Walker et al.l[l99l iVelazauez fc Whitel ll999l: iFo 



200lHKazantzidis et al.ll2007[ ). 



Font et all 



We consider the idealized case of an isolated interac- 
tion, in contrast to much work on di sk galaxy forma- 
tion done in a fu ll cosmologica l con text (iFont et aT]|2001l: 
Ardi et al.ll2003l:lBrook et al.ll2004 l200,llGauthier et al.l 



20061 [Kazantzidis et al.l 120071 )! The cosmological ap- 



proach has the relative advantage of a more realistic ac- 
cretion history. However, the isolated interactions ana- 
lyzed here offer the alternative benefit of examining the 
individual effects of a single encounter, and allow us to ef- 
ficiently sample the parameter space of interactions. Fur- 
thermore, our approach allows the simulation to be per- 
formed at much higher resolution. This helps capture the 
dynamical effects of minor mergers on the surface density 
profile at large radius, where resolution is critical. 

The simulations pre sented in this study were per- 
formed with GADGET2 (ISDringe]||2OO50 . an N-Body/SPH 
(Smooth Particle Hydrodynamics ) code using the en- 
tropy conserving formalism of ISpringel fc HernquistI 
(|2002|) . We include the effects of radiative cooling and 
star formation, tuned to fit the observed Schmidt Law 



TABLE 1 
Model Galaxy Parameters 





Sb 


Im 


M (10"' /i-i Mq ) 


95 


12 


Vc (km s-i ) 


160 


80 


c 


9 


12 


h-D {h~^ kpc ) 


4.1 


1.6 


Md/M 


0.05 


0.05 


Nh 


1.0 X 10^ 


1.3 X 10^ 


Nb 


8 X 10^ 


4.3 X 10* 



TABLE 2 
Simulation Orbital Parameters 



Name 


i 


fa 


Rp 








(ft-i kpc ) 


SbOIm30Rp1 


30° 


0.0 


5.0 


Sb2Im30Rp1 


30° 


0.2 


5.0 


Sb4Im30Rp1 


30° 


0.4 


5.0 


Sb8Im30Rp1 


30° 


0.8 


5.0 


Sb2ImORp1 


0° 


0.2 


5.0 


Sb2Im90Rp1 


90° 


0.2 


5.0 


Sb2Im150Rp1 


150° 


0.2 


5.0 


Sb2Im180Rp1 


180° 


0.2 


5.0 


Sb2Im30Rp2 


30° 


0.2 


2.5 


Sb2Im30Rp3 


30° 


0.2 


10.0 



([Schmidtl [19591: iKennicuttI I1998D . We also incorporate 
a sub-resolution multi-phase feedback model of the in- 
terstellar medium (ISM) (Springcl fc Hernquii3 llOOl) - 
softened {qEOS — 0.25) such that the mass-weighted ISM 
temperature is ^ 10^'^ - and sink particles representing 
supermassive black holes that can accrete gas and release 
isotropic thermal energy to the surrounding medium 
(Springel et al. 2005b). For furt her details on the pro - 
genitor galaxy models, we refer to lSpringel et al.l (|2005a[ ). 
and to other work done as part of a larger study of the 
effects of galaxy i nteractions on the formation and evo- 
Jiition of galaxies (Di Matteo et al.' '2005': Hopkins et alj 
2005a,b,c,d, 2006a.b,c, 2007c; Robertson ct al. 2006a,H^ 
Coxet al.ll2006allbl ld). 



A summary of the galaxy models used in the simula- 
tions is provided in Table [H including the total (baryons 
and dark matter) mass M , circular rotation velocity Vc, 
concentration parameter c, initial disk scale length ho, 
disk (stars and gas) mass fraction Md/M, number of 
dark matter particles in the halo Nh , and baryonic (stars 
and gas) particles Nb- Both are designed to be repre- 
s entative of their eponymous local bubble types (see e.g., 
Roberts fc Haynes 1994). 

The different encounter configurations considered are 
summarized in Table [H We assume zero-energy 
parabolic or bits (e = 1), as motivated by cosmologica l 
simulations ([Bensonl 120051 ; iKhochfa r fc BurkertI l2006( l. 
with radius of pericenter Rp, orbital inclination (i), and 
primary disk gas fraction fg. For the primary disk 
(Sb), we consider 20% gas disks {fg = 0.2), which 
are intended to be representative of disks in the lo- 
cal universe (jMcGaugh fc de Blok|[l997l: iBell fc de 
l2000( l. and higher gas fractions {fg = 0.4,0.8) which 
are consistent with both the more gas-rich local systems 
(iMcGaugh fc de Blokl Il99l and high redshift {z ~ 2) 
spirals ( Erb et al.ll2006f ). The secondary disk has a fixed 
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Fig. 1. — Stellar surface density profiles (S) as a function of 
radial distance from the stellar center of mass (R) for an isolated 
20% gas disk Sb2, and the remnant Sb2Im30Rp1. Included are the 
initial Sb disk (dotted), an Sb2 disk evolved in isolation for 2 
Gyr (dashed), and the Sb2Im30Rp1 remnant (solid). Also shown 
is contribution of stars from the secondary disk to the total stellar 
surface density in the remnant (dash-dot). 

gas fraction of 40% (/g — 0.4) - with the exception of 
the purely coUisionless interaction SbOImSORpI which 
has fg — - and is consistent with observations of 
dwarf galaxies and low-mass disks in the local universe 
(jSchombc rt et al.ll200ll : iGeha et ani2006( ). The initial 
spins of the two disks are not aligned, and are the same 
in all our simulations. 

3. DYNAMICAL RESPONSE OF THE STELLAR DISK TO A 
MINOR MERGER 

We find that minor mergers can create antitruncated 
stellar disks in face-on spirals. To examine this effect 
in detail and illustrate some of the generic features of a 
minor merger, we concentrate on Sb2Im30Rp1 and its 
coUisionless counterpart SbOIm3GRp1. In all cases, the 
stellar mass surface density profiles of the remnant are 
measured 1 Gyr after the final coalescence - or sev- 
eral orbital periods at the half mass radius - to allow the 
remnant disk to reach a state of approximate dynamical 
equilibrium. 

In our experiments, antitruncations are produced only 
when there is a significant supply of gas in the primary 
disk. The driving physical mechanism for producing this 
outcome represents a competition between merger driven 
inflows and transfer of angular momentum to large ra- 
dius in the remnant stellar disk; i.e., gas moves inwards 
while stars move outwards. This effect is manifest in the 
changes induced in both the gravitational potential and 
angular momentum profile during the encounter. 

In Figure [1] we present the stellar mass surface density 
profiles for Sb2Im30Rp1, compared to both the initial 
primary disk Sb2 and the same initial disk evolved in 
isolation for 2 Gyr , as a function of radial distance 
from the stellar center of mass R. The primary disk is 
stable; when it is evolved in isolation over several orbital 
periods, the inner scale length is only marginally shorter 
owing to preferential star formation occurring near the 
center (SFR - pi *: lKennicuttl[l99l . At large radius, 
there is some fluctuation of the evolved, isolated disk 
about the initial stellar mass surface density. This owes 
to Poisson noise arising from low particle counts at large 
R and the development of spiral structure from numerical 
noise associated w ith the discretized dark matter halo 
(|Hernquistlll993al) . 




10 20 30 40 10 20 30 40 
R (h"' kpc) R (h"' kpc) 



Fig. 2. — Same as Figure [T] but split into the two component 
stellar populations: "progenitor" stars which are present at the 
start of the simulation, and "new" stars which are formed during 
the interaction. The dotted line in the right panel shows the initial 
gas distribution. 

The surface density profile of the Sb2Im30Rp1 merger 
remnant shown in Figure [1] displays three key features. 
First, within 4 /i"^ kpc , the surface density profile is 
steep, indicative of a bulge component produced by the 
merger. Second, the surface density profile from 4 — 20 
kpc is nearly identical to the primary disk. Third, 
beyond 20 kpc , there is a clear excess of surface 
density relative to the initial disk of the primary. 

The excess surface density in the outer profile of the 
remnant Sb2Im30Rp1, i.e., the antitruncation of its 
disk, qualitatively - and, as we will see in § [5] and § [6l 
quantitatively - resemb les the broken exponentials of 
iPohlen fc Truiillol (|2006D . with an inner scale length close 
to that of the initial disk. In Figure [2l we separate out 
the "progenitor" stellar particles that are initialized with 
the disks and "new" stellar particles formed from the gas 
during the interaction. We find that the antitruncation 
is dominated by progenitor stars. An antitruncated disk 
in Sb2Im30Rp1 is produced independent of a fitted pro- 
file; a robust result with respect to any fitting procedure. 
Separating out the stellar particles that originate in the 
secondary stellar disk (see Figure [T]), we find that the in- 
creased outer surface density {R >2Qh~^ kpc ) is dom- 
inated by progenitor stars from the primary disk that 
have been transferred to larger radius by the interaction. 
Furthermore, we find that these large R features in the 
profile are rotationally supported - their median circular 
velocity in circular annuli is ~ 0.8 — 0.9 times the Keple- 
rian orbital velocity at that radius - and thus long-lived. 
We confirm this by evolving our remnant in isolation for 
~ 10^" years, and find that the antitruncation is not a 
transient feature. 

The surface density profiles shown in Figure [3] demon- 
strate that the coUisionless interaction SbOIm30Rp1 dis- 
plays no antitruncation in the stellar mass surface den- 
sity profile of its remnant. Rather, its surface density 
profile, which has been tilted, increasing the scale length 
at all radii, does not have a well-defined break. Features 
at large R are, as in the previous case, also rotationally 
supported. 

Although there is some evidence for bulge formation 
in Sb2Im30Rp1 at i? 2 h-^ kpc , SbOIm30Rp1 shows 
a much more pronounced bulge which remains promi- 
nent out to larger radii. This is expected; phase space 
conservation in a coUisionless interaction leads to lower 
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Isolated {t=2) 
Remnant (i = 30 , f, = 0.0) 



10 20 30 40 

R (h"' kpc) 

Fig. 3. — Stellar surface density profiles (S) as a function of 
radial distance from the stellar center of mass (R) for an isolated 
coUisionless disk SbO, and the coUisionless remnant SbOIm30Rp1. 
Included are the initial Sb disk (dotted), an SbO disk evolved in 
isolation for 2 Gyr (dashed), and the SbOImSORpI remnant 

(solid). The grey dashed line is an exponential fit to the collision- 
less remnant SbOIjviSORpI for ij > 10 kpc , meant to highlight 
that this profile, in contrast to the 20% gas remnant Sb2Im30Rp1 
(see Figure [TJ , would not likely be observed as significantly an- 
titruncated. 
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Remnant (i = 30 , f. = 0.2) 
Remnant (i = 30 , f„=0.0) 



10 20 30 40 

R (ti"' kpc) 

Fig. 4. — Gravitational potential relative to the initial disk, as a 
function of radial distance from the stellar center of mass (ij), for 
the Sb2Im30Rp1 remnant (sohd) and the coUisionless SbOIm30Rp1 
remnant (dashed). We find that dissipation efficiently deepens the 
central potential of the remnant relative to the initial disk. 



phase space densities in the core of the coUisionless rem- 
nant, and accordingly a more diffuse spheroid component 
([Hcrnquist ct al. 1993). While bulge growth via minor 
mergers is a topic worthy of further study, we postpone 
a detailed analysis to future work and point out the qual- 
itative difference between the surface density profiles of 
Sb2Im30Rp1 and SbOImSORpI. 

To investigate the physical processes driving the an- 
titruncation in our simulations, we first consider the 
gravitational potential as a function of radial distance 
from the stellar center of mass of SbOImSORpI, and 
Sb2ImS0Rp1 as compared to the initial disk in Fig- 
ure m For certain orbits, a minor merger can drive nu- 
clear inflows of gas (Hcrnquist 1 9891 : iMihos &: Hernguls^ 
Il994t iHernquist fc Mihos 1995). fuehng a centrally con- 
centrated starburst and creating a deeper potential well 
there owing to the effects of gas disspation (see Figure[2]). 
This deeper potential will contract the remnant profile, 
counteracting the broadening of the profile owing to an- 
gular momentum transfer and maintaining an inner scale 
length similar to the initial primary disk. Furthermore, 
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Fig. 5. — Total angular momentum in the stellar disk, as a func- 
tion of radial distance from the stellar center of mass (i?). Shown 
are: the initial disk (dotted), an isolated disk evolved for 2 
Gyr (dashed), and (solid) the remnant for both a gas fraction of 
/g = 0.2 (left: Sb2 disk and Sb2Im30Rp1 remnant) and a coUision- 
less interaction (right: SbO disk and SbOIm30Rp1 remnant). We 
find that angular momentum is transferred to the outer disk, which 
increases its scale length relative to the inner disk, and creates the 
antitruncation. 

the newly formed stars will be more concentrated than 
the progenitor stars (see Figure [2]), which will also tend 
to contract the inner scale length and populate the stellar 
mass surface density at small R. 

At the same time, the interaction transfers angular 
momentum and stellar mass to the outer disk. In Fig- 
ure O we show the total angular momentum in circu- 
lar annuli as a function of R. In both Sb2ImS0Rp1 
and SbOImSORpI, the angular momentum at large R 
is nearly double that of the initial and evolved disks; 
this transfer occurs whether or not gas is included. In 
Sb2ImS0Rp1, the magnitude of the angular momentum 
in shells at i? < 10 h^^ kpc is also higher, owing to the 
more efficient infiows generated by dissipation during the 
interaction. In the coUisionless case, because the inner 
potential is not as deep, the disk expands more uniformly 
in response to this transfer, and therefore does not show 
an antitruncation. Thus, the antitruncation at large R 
results from e xpansion of the oute r disk - similar to that 
first noted by 'Quin n et al.l (|1993f l - in response to a net 
transfer of angular momentum. When gas is present, the 
inner potential is deep enough to contract this inner pro- 
file and maintain an inner scale length similar to that of 
the initial primary disk. 

4. FITTING THE SURFACE DENSITY PROFILE 

In § [Sl we find antitruncated stellar disks independent 
of the fitted surface density profile. However, to facilitate 
comparison to the observational constraints on antitrun- 
cated disks, we fit profiles to the stellar mass density 
profiles of our remnants, projected face-on, and include 
both progenitor and new stars. The fit is performed, as in 
§[3l 1 Gyr after the final coalescence, so the remnant 
reaches a st ate of approximate dynam ical equilibrium. 

Following [Pohlen fc Trujillol (|2006D . we mask out the 
inner 5 kpc of the dis k and first fit a n exponen- 
tial profile (|Pattersonl Il940t iFreemanI Il97tf) with scale 
length hoi to the inner disk, then fit a second exponen- 
tial profile with scale length h 02 to the outer disk^ . The 

^ The fitting range for the outer disk was set by an initial guess 
for Rfj^. We found that the fitted value for R^,^ and the outer disk 
parameters were largely insensitive to this choice. 
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TABLE 3 
Fitted Surface Density Profiles 



Name 




h]j2/ ^Dl 






(h^^ kpc ) 


(h-^ kpc ) 




SbOIm30Rp1 


5.99 






Sb2Im30Rp1 


4.00 


1.69 


3.27 


Sb4Im30Rp1 


3.54 


1.85 


3.33 


Sb8Im30Rp1 


3.95 


1.26 


3.85 


Sb2ImORp1 


3.38 


1.72 


2.47 


Sb2Im90Rp1 


4.10 


1.24 


2.64 


Sb2Im150Rp1 


4.64 


1.61 


5.45 


Sb2Im180Rp1 


4.22 


1.27 


3.49 


Sb2Im30Rp2 


3.70 


1.80 


3.10 


Sb2Im30Rp3 


4.28 


1.61 


2.87 



break radius Rtr is defined at the intersection of the in- 
ner and outer disk profiles. The fits are performed using 
a Levenberg-Marquardt least-squares minimization rou- 
tine, with bins weighted by the Poisson error (tXi ~ E^). 
The results are tabulated in Table [31 We note as a caveat 
that the values of the parameters in our fits are somewhat 
sensitive to the manner in which the bins are weighted. 
However, using slightly different weights, such as "flux" 
(Sr^), does not qualitatively affect our results. 

Since we are considering the stellar mass distribution 
out to large radius, resolution effects are particularly 
important. We note that all of our simulations have 
more than 1.5 x 10^ stellar particles in the "outer" disk 
{ihui < R < lOhoi), and a majority have ^ 3.0 x 10^ 
over the same range. Furthermore, tripling the number 
of particles did not change the surface density profile at 
large radius {R <; 20 h^^ kpc ) by more than 15%. There- 
fore, we find that our resolution is sufficient to make ro- 
bust claims about the the stellar surface mass density 
profiles at large radius. 

5. DEPENDENCE ON THE GAS CONTENT 

Because antitruncation appears to be a dissipational 
effect, we expect the degree and location of the break to 
depend on the gas content of the primary disk. There- 
fore, we perform a set of experiments varying the gas 
fraction of the primary disk, while holding the orbital 
parameters fixed. In Figure [SI we present the stellar sur- 
face mass density profiles of the remnant, including both 
the progenitor and new stellar particles, for the four dif- 
ferent gas fractions listed in § [2] In addition, we show 
the broken exponential disk profiles listed in Table [31 

We find that SbOIm30Rp1 - the collisionless interac- 
tion, see also F igure [31- is not wel l-fitt ed by a broken ex- 
ponen tial as in lErwin eTall (120051) and lPohlen fc Truiillol 
(|200(l), and therefore would not be observed to be an- 
titruncated. Rather, the scale length increases to hoi ~ 
6h~^ kpc , relative to ho = AAA kpc initially, with 
a substantial bulge-l ike component following an i?^/'' 
()de Vaucouleursl[l959[ ) profile, which is shown in the fit 
presented in Figure [6l 

Figure [3 shows that the inner potential is deeper - and 
therefore the inner scale length shorter - for increased gas 
fractions (see Table [3]). The antitruncation is strongest 
at fg = 0.2 and 0.4, while for the highest gas fraction 
fg = 0.8, the relative scale length hiji/ho2 is signifi- 
cantly flatter, resulting in a less pronounced break. This 
is related to the composition of the initial disk; very gas 



rich disks have fewer stars remaining in the outer disk 
after the initial inflow of gas. Therefore, as the outer 
disk expands, less stellar mass resides at large radii. 

6. DEPENDENCE ON THE ORBITAL PARAMETERS 

Since the processes that create the antitruncation are 
dynamical, the degree and location of the break should 
be sensitive to the orbital parameters of the encounter. 
To investigate this, we vary the orbital parameters - the 
orbital inclination i and radius of pericenter Rp - holding 
the gas content of the primary disk fixed. 

6.1. Orbital Angular Momentum 

To test the ability of the deeper potential to mitigate 
against the expansion of the inner disk owing to angular 
momentum transfer, we vary the total angular momen- 
tum of the secondary's orbit, while fixing the orbital in- 
clination. This is done by adjusting the location of the 
pericenter of the secondary's orbit. 

We find that the antitruncation is largely insensitive 
to increasing the orbital angular momentum. Figure [8] 
shows the surface density profiles for three different radii 
of pericenter {Rp) spanning a factor of four in orbital an- 
gular momentum. The broken exponential fits are listed 
in Table [31 Though the inner scale length does increase 
over this range - despite a deeper inner potential (see 
Figure [5]) ~ it does so by less than 20%. Over the same 
range, the break radius and relative scale lengths de- 
crease by less than 10%. Therefore, we expect that at 
fixed inclination, most orbits would create similar an- 
titruncations. 

6.2. Orbital Inclination 

Varying the orbital inclination of the interaction in- 
troduces two competing effects. First, prograde minor 
mergers are more efficient than retrograde mergers at 
coupling to the rotation of the primary disk. At the 
same time, coplanar minor mergers are more efficient at 
transferring angular momentum to the stellar orbits, in- 
ducing bar formation, and centrally concentrating gas 
and stars. As a result, the closer the interaction is to 
coplanar, the deeper the remnant's inner potential. We 
present results for five different inclinations, as outlined 
in §[11 and Table H in Figures [TOl and [IH 

We first consider the prograde and polar interactions - 
SbOImORpI, SbOImSORpI, and SbOIm90Rp1 - to il- 
lustrate the combined effects of the merger-driven in- 
flow. As the inclination increases from i = Q (copla- 
nar prograde) to z = 90 (polar), the potential at small 
scale radius is shallower (see Figure [11]) . Accordingly, 
the inner scale radu of SbOImORpI, SbOImSORpI, and 
SbOIm90Rp1 are successively larger with flatter relative 
outer to inner scale lengths hD2/hDi- As a result, the 
antitruncation is less pronounced for increasing i. 

When the interaction is retrograde, as with 
Sb150ImORp1 and SbISOImORpI, the orbital fre- 
quency of the secondary is out of resonance with the 
orbits of particles in the inner disk. This makes the 
encounter much less efficient at transferring angular 
momentum to the outer disk. As a result, though 
there are antitruncations in our retrograde experiments, 
they are far less pronounced than those in prograde 
encounters. Therefore, our simulations suggest that pro- 
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Fig. 6. — Stellar mass surface density of the remnant as a function of radial distance from the center of mass R in units of the inner scale 
length hr>i (see Table [3l for fitted values), and its dependence on the initial gas fraction of the primary disk fg: from left to right fg = 0.0 
(SbOIm30Rp1), 0.2 (Sb2Im30Rp1), 0.4 (Sb4Im30Rp1), and 0.8 (SbSImSORpI). We include both the binned simulation data (crosses), and 
fitted disk profiles (solid line ). The coUisionless (fg = 0.0) remnant is not well— described by a double exponential, and instead has been fit 
with a combination disk and lde VaucouleursI II1959I ) profile. 

flows, antitruncated stellar disks are produced only when 
gas dissipation and star formation are included. More- 
over, the magnitude and location of antitruncation is re- 
lated to both the gas content of the primary disk and 
the orbital parameters of the interaction. These features 
in the surface stellar mass density profile are rotation- 
ally supported, and therefore long-lived and likely to be 
observed in local spirals. 

This merger-driven scenario for the production of an- 
titruncated disks is supported by observations of face- 
on spirals that find antitruncated disks occur more fre- 
quently in e arlier-type spirals and i n higher density en- 
vironments (jPohlen fc Truiilloll200^ . The authors note 
that the frequency of antitruncated disks increases from 
20% in Sd types to 50% in Sb types, while the fraction 
of classically truncated disks decreases from 40% to 10% 
over the same range. This agrees qualitatively with mi- 
nor mergers as the physical mechanism driving disk an- 
titruncation: spirals in higher density environments are 
more likely to have undergone minor m ergers which cre- 
ate sy stematically earlier Hubbl e tvpes (iNaab fc: BurkertI 
20031). At the same tim e, both lErwin et al.l ()2005f ) and 




10 20 30 40 

R (h"' kpc) 

Fig. 7. — Response of the gravitational potential to a mi- 
nor merger as a function of radial distance from the stellar cen- 
ter of mass R, and its dependence on the initial gas fraction of 
the primary disk fg-. fg = 0.0 (dashed, SbOImSORpI), 0.2 (solid, 
Sb2Im30Rp1), 0.4 (dash-dot, Sb4Im30Rp1), and 0.8 (long dash, 
Sb8Im30Rp1). The dotted line shows the initial potential of the 
primary disk. 

grade minor mergers will be most effective at producing 
antitruncations. 

7. DISCUSSION 

We find that minor mergers can create antitruncated 
stellar disks in face-on spiral galaxies, and that this 
antitruncation is produced by a competition between 
merger driven inflows of gas and transfer of angular mo- 
mentum to large R in the remnant stellar disk that moves 
stars outwards. Because this process requires gaseous in- 



Pohlen fc Truiillol ()2006[ ) show observational evidence for 



asymmetries or recent interactions in antitruncated sys- 
tems, which further supports a merger-driven scenario. 

Using the extended Press-Schecht er formalism 
( Jenkins et al .1120011 ) and the method of iLacev fc Cold 
( 1993f ) to estimate hal o merger h i storie s, and assuming 
the cosmology of Soe rgel et al.l ()2003[ ). we find that 
SB-type halos {Mto't 10^^ h'^ Mq ) are likely to 
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Fig. 8. — Stellar surface mass density of the remnant as a function of radial distance from the center of mass R in units of the inner scale 
length hui (see Table[3]for fitted values), and its dependence on the radius at pericenter Rp of the orbit of the secondary: from left to right 
Rp = 2.5 {Sb2Im30Rp2), 5.0 (Sb2Im30Rp1), and 10.0 (Sb2Im30Rp3) h^^ kpc . Wc include both the binned simulation data (crosses), and 
fitted disk profiles (solid line). 

IPohlen fc TrlluIIol ((200l) and lErwin et all ((2005h . We 

find, however, that the break radius in our simulations is 
at the low end of the observed range : IPohlen fc Truiillol 
(|2006f l and lErwin et all (|2005D find Rbr/hoi = 3-6 
while in our simulations Rbr/hoi = 3 — 4. This may 
owe either to: (1) the s malle r mass of a typical galaxy in 
the IPohlen fc Truiillol ()2006D sample, or (2) the hmited 
range of parameter space spanned by our s imulations. 
The observations of IPohlen fc Truiillol (|2006D appear to 
be dominated by somewhat lower-mass spirals than the 
Milky- Way mass primary disk in the interactions exam- 
ined here. This c ould p otentially lead to shorter inner 
scale lengths (e.g.. ICour tcau 1996; dc Jong 1996), which 
may tend to increase the average observed Rbr/hoi ra- 
tio. Also, our simulations sample only a small subset 
of the parameter space for individual interactions. The 
real merging history of galaxies likely involves a variety 
of mass fractions and multiple mergers that may produce 
subtly different effects. However, we find that our simu- 
lations of minor mergers generically lead to antitruncated 
disks for a range of orbital geometries, and therefore rep- 
resent a viable mechanism for producing these features. 

More locally, Ibata et all ([2005) and Ibata et al. (200'^ 
recently observed an extended stellar disk in M31. 
Though it is not entirely clear whether or not this fea- 
ture is preceded by a well defined break in the sur- 
face brightness profile, it is possible that the extended 
dis k represents a n ant itrunca -tion of the type observe d 
bv lErwin et al.l ^2005h and IPohlen fc Tr^aiilloj HOO^. 
Ilbata et al.l ( 20051 ) estimate that it contains roughly 10% 



10 20 30 40 

R (h"' kpc) 

Fig. 9. — Response of the gravitational potential to a minor 
merger as a function of radial distance from the center of mass 
R, and its dependence on the radius at pericenter Rp of the orbit 
of the secondary: Rp = 2.5 (dashed, Sb2Im30Rp2), 5.0 (soUd, 
Sb2Im30Rp1), and 10.0 (dot-dash, Sb2Im30Rp3) h'^ kpc . 

experience of order one 1:8 minor merger from z — \ 
to the present-day (| Hopkins et al.l l2007ai rb[). In our 
simulations, strong antitruncations are produced when 
the orbit of the secondary is inclined (0° ^ i ^ 90°), pro- 
grade, and has moderate angular momentum {Rp hu). 
So, if all orbits are distributed isotropically - i.e., equally 
likely in bins of dcosi - and follow the distribution of 
Rp a s inferred fro m cosmological N-body simulations 
(e.g.. lBensoi3l2005D . then we would expect ^ 40% - 50% 
of Sb type spirals to have pronounced antitruncations. 

At the same time, our fits agree quantitatively 
with the relative scale lengths /i_D2//ini observed by 
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Fig. 10. — Stellar surface mass density of the remnant as a function of radial distance from the center of mass R in units of the inner 
scale length hui (see Table[3]for fitted values), and its dependence on the orbital inclination i of the interaction: from left to right i = 
(coplanar prograde, Sb2ImORp1), 30 (prograde, Sb2Im30Rp1), 90 (polar, Sb2Im90Rp1), 150 (retrograde, Sb2Im150Rp1), and 180 (coplanar 
retrograde, Sb2Im180Rp1). We include both the binned simulation data (crosses), and fitted disk profiles (solid line). 
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Fig. 11. — Response of the gravitational potential to a minor 
merger as a function of radial distance from the center of mass R, 
and its dependence on the orbital inclination i of the interaction: 
1 = (coplanar prograde, dashed, Sb2ImORp1), 30 (prograde, solid, 
Sb2Im30Rp1), 90 (polar, dot-dash, Sb2Im90Rp1), 150 (retrograde, 
double dot-dash, Sb2Im150Rp1), and 180 (coplanar retrograde, 
long dash, Sb2Im180Rp1). The dotted line represents the initial 
potential of the primary disk. 

of the stellar mass and 30% of the angular momentum of 
the total stellar disk, as compared to 5% of the mass and 
45% of the angular momentum in the antitruncated disk 
of Sb2Im30Rp1. Furthermore, the kinematics - specif- 
ically, the dispersion of circular velocity lags relative to 
Keplerian rotation - of resolved stars in the extended disk 
show evidence of dynamical heating which could have 
been caused by a minor merger (e.g. Quinn " eiTanflOQl 
IWalker et all Il996t IVelazguez fc White ,19991 1 There- 
fore, although we cannot say with certainty that the 



extended disk of M31 represents a local example of an 
antitruncation, it is broadly consistent with our model- 
ing. 

8. CONCLUSION 

We use hydrodynamic simulations to investigate minor 
mergers as a physical mechanism for creating antitrun- 
cated disks in face-on spirals. We find that the antitrun- 
cation is produced by two competing effects: merger- 
driven gas inflows deepen the central potential and con- 
tract the inner profile, while at the same time angular 
momentum is transferred to large radius and causes the 
outer disk to expand. Because the inflows are far more 
efficient when gas dissipation is included, the antitrun- 
cation is produced in our experiments only when a sig- 
nificant gas supply is present in the initial primary disk. 
This effect is also only seen when the interaction is pro- 
grade, rather than polar or retrograde, with moderate 
{Rp ~ ho) orbital angular momentum. 

Our merger-driven scenario for producing antitrun- 
cated disks yields results that agree with obser- 
vatio ns of local face- on spirals (jErwin "eTaLl [20051 : 
Pohl en fc Truiillol l200(if ). both in terms of the parame- 
ters of the antitruncation and its frequency with Hubble 
type. Therefore, we find that minor mergers are a viable 
physical mechanism for producing antitruncated disks. 
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